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Abstract Following the reports of candidate chiral dou- 
blet bands observed in cesium isotopes, the possible chiral 
candidates and the evolution of three-dimensional rotation 
in 1207154Cs are investigated within the microscopic three- 
dimensional tilted axis cranking covariant density functional 
theory (3DTAC-CDFT). By investigating the evolution of 
the polar angle 0 and azimuth angle @ as a function of rota- 
tional frequency fa, the transition from the planar rotation 
to the chiral rotation has been found in !?! ^5? Cs, The corre- 
sponding critical rotational frequency ox; of the appearance 
of chiral aplanar rotation decreases as neutron number in- 
creases, which can be attributed to the neutrons in (gd) and 
(sd) shells having smaller angular momentum components 
along both the short and long axes, and larger components 
along medium axis, respectively. In comparison, only pla- 
nar rotation has been obtained in !?®!34Cs. With these inter- 
pretations, the obtained J ~ œ and energy spectra as well 
as B(M1)/B(E2) values show reasonable agreement with 
the available experimental data. In addition, the evolution 
of quadrupole deformation f and triaxial deformation y are 
also discussed. 


1 Introduction 


Chirality has attracted general interest in natural sciences, 
such as chemistry, biology, and physics. In nuclear physics, 
chirality, i.e., the spontaneous chiral symmetry breaking, 
was originally suggested by Frauendorf and Meng in 
1997 [1]. In a triaxially deformed nucleus with the unpaired 
particle(s) and hole(s) in the high-j orbit, the particle(s) 
and hole(s) align their angular momentum vectors along the 
short(s) and long(l) axes, respectively, while the collective 
core aligns along the medium(m) axis, the total angular mo- 
mentum vector lies outside the three principal planes in the 
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intrinsic frame (aplanar rotation) with left- or right-handed 
orientations, which can form a chiral system. The restoration 
of the spontaneous chiral symmetry breaking in the labora- 
tory frame could give rise to the observation of the so-called 
chiral doublet bands, i.e., a pair of nearly degenerate AJ = 1 
bands with the same parity [1]. In fact, a pair of AJ = 1 bands 
found in '34Pr with the zh), /2 8 Why; /2 configuration ob- 
served early in 1996 [2] have been reinterpreted in Ref. [1] 
as a candidate for chiral doubling. Thereafter, in 2001 sim- 
ilar low-lying doublet bands were reported in 55Cs, 57La, 
and 6;Pm N = 75 isotones of !?^Pr, and an island of chiral 
rotation was suggested in the A ~ 130 mass regionin [3]. It 
should be noted that clearly identifing the chiral motion is 
difficult, for example the chirality in !?^Pr aforementioned 
was later rejected in 2006 [4, 5], and in 2011 [6] other chiral 
candidates in '?^Pr were proposed. 

Since this pioneering work, a lot of experimental and 
theoretical works related to nuclear chirality have been done, 
see Refs. [7-10], for a brief review. Thus nuclear chiral- 
ity becomes a hot topic and related issues such as MyD 
(multiple chiral doublet bands in one nucleus) [11—15], the 
nuclear Chirality-Parity (ChP) violation [16, 17] and so on 
have been widely discussed. 

Up to now, more than 60 candidate chiral nuclei have 
been reported in the A ~ 80,100,130,160,and 190 mass 
regions of the nuclear chart, see, e.g., data tables [18]. In 
particular, the reported candidate chiral nuclei in A~130 
mass region form a large chiral island, where the ce- 
sium isotopes have the most chiral candidates. So far, chi- 
ral doublet bands have been observed in the odd-odd nu- 
clei 122-124,126,128,130,132¢, [3, 19-24] with configuration 
Thy, NEVA» and corresponding chirality are also well 
explained by different theoretical models [19, 20, 23-38]. In 
fact, chiral doublet bands were first identified in the A ~ 130 
nuclei including 130Cs [3]. The chiral doublet bands found in 
128Cs was proposed as the best example to reveal the chiral 


symmetry breaking [19], and the succeeding g-factor mea- 
surements can give important information on the relative ori- 
entation of the three angular momentum vectors [39, 40]. 
Meanwhile, the lifetimes of excited states belonging to the 
chiral partner bands have been reported in Cs, which is 
the first time a large set of experimental transition proba- 
bilities are in qualitative agreement with all selection rules 
predicted for the strong chiral symmetry breaking limit [21 ]. 
Therefore, it is naturally interesting to find nuclear chiral- 
ity in other cesium isotopes. In neighbouring odd-odd iso- 
topes !7°:!34Cs, only one rotational band with the config- 


uration zh, phy» was reported [32, 41]. In the heav- 


ier '34Cs, the experimental observations [32, 42] indicate 
that the structure of the band is distinctly different from its 
lighter odd-odd isotopes, and was interpreted as a possible 
magnetic rotation band [42]. Therefore, it is worthwhile to 
recheck the rotational modes in 1?®134Cs and the possible 
boundaries of the cesium isotopes. 


For the odd-A cesium isotopes, doublet bands involving 
the configuration component Thy, phy, with the third 
neutron quasi-particle or hole by the g7/2/d5/z or 8} /2/d3/2 
orbits in 145095 05,07, 09,1310 [43-48] have been observed. 
Theoretically, the constrained triaxial relativistic mean-field 
(RMF) approach [49, 50] has been performed to study the 
possible chirality in odd-A cesium isotopes !2?:127.29.131 Cs, 
and the existence of chirality is demonstrated and expected 
based on different high-j particle-hole configurations and 
triaxial deformations. In addition, a rotational band with the 
configuration component zh}, p vhi, (sd) ! has been 


observed in !??Cs [51]. Therefore, it is also interesting to 


investigate the chirality in other odd-A cesium isotopes, i.e., 
121,123,133. 


In order to describe the chiral geometry microscopi- 
cally, we adopt the three-dimensional tilted axis cranking 
based on the covariant density functional theory (3DTAC- 
CDFT). This approach has been developed [52] and used to 
investigate the nuclear chirality in 106Rh [52], 106 Ag [53], 
135Nd [54] and !°2-!07Rh [55]. In these studies, the exis- 
tence of the critical frequency @cir corresponding to the 
transition from planar to aplanar rotation is supported. For 
example, in Ref. [54], a classical Routhian was extracted by 
modeling the motion of the nucleons in rotating mean field 
as the interplay between the single-particle motions of sev- 
eral valence particle(s) and hole(s) and the collective motion 
of a core-like part. This classical Routhian gives qualitative 
agreement with the 3DTAC-CDFT result of @,;,;. In this pa- 
per, following the similar procedures outlined in Refs. [52— 
55], we will pay attention to chirality in 197 1?^Cs within 
3DTAC-CDFT. 


2 Theoretical framework 


The 3DTAC-CDFT based on the zero-range point-coupling 
interaction has been successfully applied to describe the chi- 
ral rotation [52-55]. The detailed formalism of the 3DTAC- 
CDFT can be found in Ref. [52]. Here, a brief introduction 
is presented. 

The CDFT starts from a Lagrangian, and the correspond- 
ing Kohn-Sham equations have the form of a Dirac equa- 
tion with effective fields S and V derived from this La- 
grangian [56-60]. In the 3DTAC method, these fields are 
calculated in the intrinsic frame rotating with a constant an- 
gular velocity vector @, pointing in an arbitrary direction in 
space [52] 


[a -(p—V) 4- B(m- S) - V — o-J]w. = Eryr. (1) 


Here, J is the total angular momentum of the nucleon 
spinors, and the relativistic fields S, V, and V read 


S(r) = asps + Bsps + Ysps + s^ ps, 

V°(r) = av py + wey + ôv APv + etry pry 
+36ryApry +eA°, (2) 
V(r) = ay jy +W Gv) + vA Sy + BOrViry 

+3 6rvAjry +eA. 


The S, V, and V are the relativistic scalar field, the time- 
like component of vector field, and the space-like compo- 
nents of vector field, respectively, which are in turn cou- 
pled with the nucleon densities and current distributions. In 
the above equation, the ps, py, Prv, jv. jrv represent var- 
ious densities and currents, and the as, Bs, Ys, Os, Ov, Ww, 
ôv, Ory, Ory are nine parameters of the Lagrangian. The 
Dirac equation is solved in a set of three-dimensional har- 
monic oscillator basis iteratively, and one finally obtains the 
single-nucleon spinors Wg, the single-particle Routhians €x, 
the total energies, the expectation values of the angular mo- 
menta, transition probabilities, and so on. Using the semi- 
classical cranking condition (J) - (J) = I(I + 1), one can re- 
late the magnitude of the angular velocity @ to the angular 
momentum quantum number J. Meanwhile, the orientation 
of @ is determined self-consistently by minimizing the total 
Routhian. 

From the Dirac equation, the physical observables, in- 
cluding the quadrupole moments and magnetic moments, 
and the electromagnetic transition probabilities B(M1) and 
B(E2), can be calculated. 

The quadrupole moments Q20 and Q22 are 
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Fig. 1 (color online) Single proton (a) and neutron (b) Routhians near the Fermi surface in 120-194Cs as a function of the neutron number based 
on the fixed configuration when the rotational frequency is 0.3 MeV. The single proton and neutron levels belonging to the orbits 7/5 /ds/; and 
51/2/ds/» are marked by solid black and blue lines, respectively, and the levels belonging to /1; /2 orbit with negative parity are indicated by dashed 
red lines. The red, grey and blue filled circles indicate the occupied levels in the orbits 11/2, 87/2 /d5 and s, /9/d3/. 
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The atomic quadrupole deformation parameters B and 
triaxial deformation parameters y can be obtained 


B = /a)+2a5,, Y= arctan v] À (4) 
a20 


The nuclear magnetic moment in units of the nuclear 
magneton is given by, 


2 

mc t 

u=} m far Savin ay, + ky, BEY), (5) 
k>0 


where the charge q is 1 for protons and 0 for neutrons 
in units of e, K is the nucleon anomalous gyromagnetic fac- 
tor, K, = 1.793, K, = —1.193, L and X are respectively the 
orbital angular momentum and spin [61]. In the following 
calculations, the ratio mc? /hc in Eq. 5 is taken as 1, as 
Refs. [52, 62] 

The transition probabilities B(M1) and B(E2) are calcu- 
lated in the semiclassical approximation, 


B(M1) = [tuc sin 8 + cos 0 (L1, cos @ + LL sin 9) 


+ (Lycos 9 — using)’ ) 
2 


2 
B(E2) == |Q5, sin? 0 + Zo (1 + cos? 0) cos29 
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+ (Q5, cos O sin29) z 
(6) 


where B and D are the quadrupole moments of pro- 
tons, and 0 and @ are the orientation angles of the total an- 
gular momentum J in the intrinsic frame. 

In the present work, we adopt the point-coupling La- 
grangian PC-PK1 [61] and the calculations are free of ad- 
ditional parameters. The Dirac equation Eq. 1 is solved in a 
three-dimensional Cartesian harmonic oscillator basis with 
10 major shells. The pairing correlation is neglected in the 
calculations, but one has to bear in mind that the pairing cor- 
relation might have some influences on the descriptions of 
critical frequency [37] as well as the total angular momen- 
tum and B(M1) values [63, 64]. 


3 Results and discussion 


In the present study for !2°-!*4Cs, two quasi-particle config- 
urations Thy, 8 Vli and three quasi-particle configura- 


tions nhi, vhi (gd) & nhi, y@VAj (sd)! are as- 


signed for odd-odd nuclei and the odd-A nuclei, respectively, 
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Fig. 2 (color online) The evolution of the polar angle 0 and azimuth angle q for the total angular momentum J as driven by the increasing 
rotational frequency ñw for configurations TAY 585, 8Vh] j (gd)" (n = 8— 14) in 207 26Cs, zh], 585/59 Vh sd)" (n = 1—2) in ?7- 286g 


and zl, gf, VI, (sd) (n = 1.— 6) in 97 PACs, 


similar to the earlier studies in Refs. [3, 19-24, 32, 41, 43— 
47, 49—51, 65]. In the calculations, an unpaired proton is 
always occupying the bottom of the 1; /; orbit and the oc- 
cupation of the valence neutrons in /1,; /? orbit are fixed by 
tracing the single particle levels with the increasing frequen- 
cies as shown in Fig. 1, while other nucleons are treated self- 
consistently by filling the orbits according to their energies 
from the bottom of the well. The (gd) represents the low-j 
neutron orbits 187/2 and 2d5/», and (sd) represents 2d3/2 and 
351 /2, as it is difficult to distinguish them as a strong mixing 
with each other. Note that the six or eight neutrons in /1; /2 
orbit are antialigned, and the corresponding neutron config- 
uration is vhi. In this way, the valence nucleon configu- 


rations of !2°-!26Cg within the 3DTAC-CDFT calculations 

are zh 587,8 Vhi; (gd) (n = 8 — 14), of ?77 Cs are 

zh 585,8 vhi; p(s)" (n —1—2), and of I?-U^Cs are 
1 4 9 n a 

Th, 585 49 Vhi (sd) (n — 1-6). 


In the 3DTAC-CDFT calculations, the orientation of 
the angular velocity @ with respect to the principal axis 
can be determined in a self-consistent manner either by 
minimizing the total Routhian or by requiring that @ is 
parallel to the total angular momentum J at a fixed @ 
value. Here, the polar angle 0 and the azimuth angle Q 
are used to denote the direction of the angular velocity 
@ = O(sinOcos Q,sinOsin 9,cos Q). The 0 is the angle be- 
tween the angular velocity @ and the long (/) axis, while 
@ is the angle between the projection of @ onto the short- 
medium (sm) plane and the short (s) axis. It has been noted 


that Q can be used to characterize the chirality of a rotational 
system [64]. 

To examine the possible presence of chiral geometry, the 
self-consistently obtained orientation angles 0 and q of the 
total angular momentum J in the intrinsic frame are shown 
as a function of the rotational frequency 710 in Fig. 2. The 
polar angle 0 in 97 ?^Cs has similar behavior, i.e., they in- 
crease with the rotational frequency. Nevertheless, the polar 
angle 0 is always larger than 50? in all cesium isotopes. It 
is attributed to that the angular momentum alignment along 
the s axis coming from the proton particles in the hj; /2 Or- 
bit is much larger than that along the / axis from the neu- 
tron holes in the 41/2 orbit (cf. Fig. 4). For comparison, 
the azimuth angle Q for '!~!°3Cs is zero at low rotational 
frequencies, corresponding to the planar rotation in the s/ 
plane. Above the limited rotational frequency correspond- 
ing to the so-called critical frequency @,,i of chiral rotation 
[37, 66], the values of 9 become nonzero and it results in 
the transition from planar to aplanar rotation. Note that the 
kinks appear at the curves of 0 for !2!~!33Cs with the critical 
frequency. 

The azimuth angle @ is always zero for !?®!34Cs cor- 
responding to the planar rotation in the s/ plane, even 
the rotational frequency goes up to hœ = 0.55 MeV for 
120Cs and hw = 0.45 MeV for '*4Cs. By increasing ro- 
tational frequency, the results can not be converged for 
the configuration of TR} 1984 /y@VA (gd) in 9 Cs and 


Th] 585, Vli j (sd) in 134Cs. Namely, it might mean 
that /@,,it in the present 3DTAC-CDFT calculations is larger 


than ha = 0.55 MeV for !°Cs and f@ = 0.45 MeV for 
134Cs even if it exists, which corresponds to I ~ 24 h, out 
of the spin range observed currently [32, 42, 47]. Thus, the 
present theoretical analysis does not support chirality exsit- 
ing in 29 D49.. 
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Fig. 3 (color online) Composition of the total angular momentum 
at different rotational frequencies in 3DTAC-CDFT calculations for 


20 Cs and !34Cs with the configuration zh! io Vh. 


In addition, the previous study indicates a possible MR 
nature for this configuration zh, phy» in ^Cs based 
on the decreasing trend of B(M1) [42]. Composition of the 
total angular momentum at different rotational frequencies 
in 3DTAC-CDFT calculations for P?Cs and !?^Cs with the 


configuration Thy, a8 Vh have been given in Fig. 3. In 


the present 3DTAC-CDFT calculations for the 1?%134Cs, as 
the four proton particles in the g7/2 orbital are paired, the 
proton angular momentum comes mainly from the parti- 
cles in the h11/2 orbit, which aligns along the s axis. The 
neutron hole in the 4; /2 orbit and neutrons in low- j orbits 
(gd)/(sd) give substantial contributions to the neutron an- 
gular momentum, which leads to large components both in / 
axis and s axis. Higher spin states in the band are created by 
aligning the neutron angular momentum toward the s axis, 
while the proton angular momentum keeps unchanged in the 
s axis. Thus, it is probably the planar rotation in 94s, 
The chirality in nuclei with stable triaxial deformation 
is due to the aplanar rotation formed by the valence par- 
ticle(s), valence hole(s) and collective core angular mo- 
mentum vectors [1]. Therefore, an important work is the 
transition point from planar to chiral rotation, i.e., the crit- 
ical frequency Qi. As shown in Fig. 2, the calculated 
critical frequencies @ert all decrease with the increasing 
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Fig. 4 (color online) Contributions of the valence neutrons in the 
11/2, (gd) and (sd) orbits to the total angular momenta along the s, 


m and 1 axis in 124125:126C¢. 


mass number for three different configurations in 120-126Cs, 


127,28 Cs and 1297134 Cs. To understand this interesting be- 
havior, one has to analyze the angular momentum geometry. 
In 3DTAC-CDFT calculations, the total angular momentum 
comes from the individual nucleons in a coherent superpo- 
sition manner. 


For all the present cesium isotopes, only one proton par- 
ticle occupying at the bottom of the 1; ; orbit contributes 
the angular momentum of roughly 5.5 A along the s axis. It 
should be emphasized that the angular momentum contribu- 
tion from proton part is approximately along the s axis at 
low rotational frequency. Therefore, the angular momentum 
components along / axis and m axis are mainly from neutron 
part. At this point, !**Cs, !?5Cs and !7°Cs are used as the 
demonstration to study the angular momentum contributions 
of the valence neutrons in the /1;;/?, (gd) and (sd) orbits 
along the s, m and l axis, which is shown in Fig. 4. In con- 
trast, five neutron holes locating at the top of the /;, /2 orbit 
contribute about 3.5 /i to the angular momentum along the / 
axis since four of them are antialigned. The only difference 
is the neutrons distributed over the (gd) orbits and results in 
different configurations zi, a8 Vh, (gd)" (n = 12— 14) 


for 124— 1260s 


Figure 4 shows that the angular momentum increment 
of the (gd) orbits along the s and / axis will become smaller 
from !?4Cs to ?Cs and the angular momentum increment 
along the m-axis will become larger. The similar behavior 
exists with the (gd) or (sd) orbits from !7°Cs to '74Cs. So, 
we speculate that the angular momentum increment of the 
(gd) or (sd) orbits along the s and / axis will become smaller 


with the increasing of the neutron number and the angular 
momentum increment along the m-axis will become larger. 
Itis easier for the nucleus to form the chiral rotation, and the 
corresponding critical frequency @erit will decrease. There- 
fore, it is crystal clear that the critical frequency with the 
same configuration in Fig. 2 always tends to move to the 
left, which is similar to the earlier studies in Ref. [55]. It 
should be noted that the !33Cs, the (sd) orbits neutron in 
this nucleus is transformed from the particle state character 
to the hole state behavior since a different change to others. 


The critical frequency also affects the behaviors of ex- 
perimental spin-rotational frequency /(A@) relationship. As 
shown in Fig. 5, the calculated values of spin I(ħ@) agree 
well with the experimental data for 129 !?3Cs, and overesti- 
mate the data for 124—-127,.129—131133 s. In addition, the cal- 
culated values of spin Z(ñ@) is much higher than the exper- 
imental data for !28:!32:134¢Cg_ For 124-13 (^s. the effective 
neutron and proton pairing correlation and beyond mean- 
field effects which are not considered here may play the im- 
portant roles, which has been confirmed in Ref. [67, 68]. 
Furthermore, the back-bending phenomenon in !78:!32:134Cg 
which are not considered here may also affect the present 
calculated results much higher than the experimental data. 
Obviously, the calculated and experimental rotational fre- 
quencies æ for the yrast bands in !?! -??Cs behave sim- 
ilaly with respect to the angular momentum. The bands 
show nearly a straight line at low rotational frequency. A 
kink will appear when the spin reaches the critical frequency 
@erit in 1217133 Cs, The occurrence of this kink is due to the 
transition from the planar to the chiral solutions and the an- 
gular momentum component will increase along the m axis 
causing huge changes in the angular momentum. The cor- 
responding kink of experimental I(ħ@) appears almost the 
same as the calculated A@g,it results. Further experimental 
efforts for !?Cs are encouraged to verify the present con- 
clusion on the A@erit. For !2°134Cs, previous experimental 
investigation does not support static chirality interpretation 
in the spin range observed [32, 42, 47], which agrees with 
the theoretical result of no kink appearing. 


As aforementioned, the positive-parity rotational bands 
have been observed in 1297 !5^Cs in the previous works [3, 
19-24, 43-46, 49, 50], and the comparisons between the 
3DTAC-CDFT results of the excitation energies and the 
available experimental data are shown in Fig. 6. In the 
present mean-field level, either the chiral vibrations or the 
tunneling between the left- and right-handed sectors are not 
taken into account. Therefore, only the band with lower ex- 
citation energies can be reproduced. However, it can be seen 
in Fig. 6 that the experimental excitation energies of the 
lower band can be reproduced well. To describe the partner 
band, one needs to go beyond the mean-field calculations by 
combining, for example, the methods of the random phase 


approximation [35] or the collective Hamiltonian with the 
CDFT [64, 69]. 


In Fig. 7, the available experimental data of 
B(M1)/B(E2) in '7°-14Cs [3, 19-24, 43-46, 49, 50] 
are displayed in comparison with the 3DTAC-CDFT results, 
which are derived in the semiclassical approximation 
from the magnetic and quadrupole moments. Note that the 
deformation parameters (fy) change only slightly along 
the band, so the corresponding B(E2) values are almost 
constant. However, the B(M1) values decrease smoothly 
due to the continuous variation of rotational frequency and 
the closing of the neutron and proton angular momentum 
vectors, which mainly align along the short and long 
axes, respectively. This leads to the smooth-decreasing 
tendency for the B(M1)/B(E2) ratios. Therefore, the ratios 
of B(MI)/B(E2) by 3DTAC-CDFT are different from 
the experimental data, and only the smooth-decreasing 
tendency for the B(M1)/B(E2) ratios can be given. In 
particular, the magnetic moments are derived from the 
relativistic expression of the effective current operator as in 
Ref. [57]. As shown in Fig. 7, the 3DTAC-CDFT results for 
121,123,125—134Cs show good agreements with the data, while 
the calculation overestimates the data for 129.122. 24s. The 
theoretical results show that B(M1)/B(E2) for !?!-?*Cs 
have a smooth falling behavior in the planar rotation, and 
the falling tendency slows down above the critical frequency 
Ox;i. It once again proves that the nuclei will transition 
from the planar rotation to the chiral rotation above the 
critical frequency Writ. The falling tendency shows steeper 
and steeper with increasing neutron number. Further efforts 
of including the neutron and proton pairing correlations will 
be helpful to justify the present results. 


The calculated quadrupole deformation parameters f 
and triaxial deformation parameters y are given in Fig. 8. 
It could be found that B and y are stable and there is 
only one slight change with the increasing rotational fre- 
quency in !297!34Cs, With the increasing neutron number, 
the quadrupole deformation f in cesium isotopes gradually 
decreases which is consistent with the results in Refs. [70, 
711]. It could be understood that the number of neutrons in 
!20Cs is the mediant of the magic number, so the quadrupole 
deformation B will gradually decrease. For the triaxial de- 
formation y, it increases with the number of neutrons in 
120-126Cg and then decreases with the number of neutrons 
in 127-154Cs. It could be understood that the triaxial defor- 
mation y in !2°-!?6Cs gradually increases with the num- 
ber of neutrons, the neutron holes contribution becomes 
stronger, and the deformation changes from prolate ellipsoid 
to oblate ellipsoid, which leads to the triaxial deformation y 
to grow up. Similarly, the neutron particle contribution be- 
comes weakened with the increasing number of neutrons in 
127-1905. and the deformation evolves from oblate ellipsoid 
to prolate ellipsoid, which leads the triaxial deformation y 
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Fig. 5 (color online) The spin as functions of the rotational frequencies 7p by 3DTAC-CDFT in comparisons with available experimental 
data [3, 19-24, 32, 41, 43-51] in 1207 4C, 
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Fig. 6 (color online) The calculated rotational excitation energies in !?°-!34Cs by 3DTAC-CDFT, as a function of the spin in comparison with 
the data for the chiral bands observed in Refs. [3, 19—24, 32, 41, 43—51]. 


to decline. However, the suddenly increasing triaxial defor- ^ the planar rotation of zero azimuth angle @ has been ob- 
mation y between ?Cs and ??Cs is due to the increasing tained. Therefore, the existence of chirality can not be com- 
neutron number in 71; j» orbit. pletely obtained based on the high-j particle-hole configu- 


rations and stable triaxial deformation [49]. 
It could also be found from Fig. 8 that the triaxial de- 


formation in !*4Cs is approximately zero, and it further in- 
dicates there is no chirality in '*4Cs, in agreement with the 
previous result [32]. The triaxial deformation parameter y 
in Cs show relatively obvious value (~ 15°), although 
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Fig.7 (color online) The calculated B(M1)/B(E2) ratios in !2°-!34Cs by 3DTAC-CDFT, as a function of the spin in comparison with the data 


for the chiral bands observed in Refs. [3, 19—24, 32, 41, 43-51]. 


60 


Fig. 8 (color online) The evolutions of the deformation parameters f 
and y in the 3DTAC-CDFT calculations for !20—15^Cs with the config- 


uration 7h}, BVI (gd) and mh}, n2Vhii/ (sd)", respectively. 


4 Summary 


In summary, a fully self-consistent and microscopic inves- 
tigation for the evolution of chirality in the cesium isotopes 
120—134Cs has been performed with the 3DTAC-CDFT. By 
investigating the evolution of the polar angle 0 and azimuth 
angle @ for the total angular momentum J as driven by 
the increasing rotational frequency ha, the transition from 


the planar rotation to the chiral rotation has been found in 
121-133 C$. In comparison, only planar rotation has been ob- 
tained in !?®!34Cs, and it is probably not the magnetic rota- 
tion. 


Moreover, the critical frequency Merit in 121—133Cs de. 
creases with the increasing neutron numbers, and it also can 
be observed by the behaviors of J and B(M1)/B(E2) vs. ha 
relationship. Taking !*4Cs, !?5Cs and !*°Cs as examples to 
study the angular momentum contributions of the valence 
neutrons in the 11/2, (gd) and (sd) orbits along the s, m 
and / axis, the angular momentum increment of the (gd) or 
(sd) orbits along the s and / axis will become smaller with 
the increasing of the neutron number and the angular mo- 
mentum increment along the m-axis will become larger. It 
is easier for the nucleus to form the chiral rotation, and the 
corresponding critical frequency @,,;, will decrease. 


In addition, the calculated 7 ~ h@ show good agree- 
ment with the available experimental data for 1297 !?3 Cs, 
and overestimate the data for 1^- !34Cs. The calculated en- 
ergy spectra for yrast bands of !?°-!*4Cs are reasonable 
with the corresponding experimental data. The 3DTAC- 
CDFT results show a good agreement with the data of 
B(M1)/B(E2) for 21.123.125-19 Cs. while overestimates the 
data for !20!22.124Cs, The quadrupole deformation param- 
eters B and triaxial deformation parameters y in !2°-!74Cs 
show obvious regular changes. The f) decreases in !2°-!4Cs 
with the increase of neutron numbers. y increases in 
120-2605 and decreases in !?7-!9?^Cs, which is related to 
the occupation of neutron numbers. Triaxial deformation in 
0 Cs and !34Cs indicates that the existence of chirality can 
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not be completely obtained based on the high- j particle-hole 
configurations and stable triaxial deformation. 

Recently, the chiral bands of !!°Cs [72] based on config- 
urations with only protons are observed for the first time in 
the A ~ 130 mass region, i.e., a configuration with three pro- 
tons, one in the strongly coupled [404]9/2* orbital which 
does not change orientation with increasing rotational fre- 
quency, and two in the 41/2 orbital which reorients to the 
rotation axis, and it is different from the common chiral ge- 
ometry picture from both protons and neutrons. It is an ex- 
cellent significant and interesting work. The preliminlar re- 
sults of 3DTAC-CDFT calculation don't support the exsi- 
tence of aplanar rotation with such configuration. At the 
same time, not only the TAC-CDFT, but also more other 
theoretical models such as particle-rotor model is helpful to 
justify the present results and understand the underlying ro- 
tational structure in !?Cs. Investigations in these directions 
are in progress. In addition, more efforts on the investiga- 
tion of nuclear chirality in the A ~ 130 mass region are still 
needed, including different theoretical methods, the evolu- 
tion of the isotones, possible configurations and so on. 
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